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Abstract 

The characterization of small samples of cold bosonic atoms in rotating microtraps has recently 
attracted increasing interest due to the possibility to deal with a few number of particles per site 
in optical lattices. In this paper we consider two-dimensional systems of few cold Bose atoms 
confined in a harmonic trap in the XY plane, and submitted to strong rotation around the Z 
axis. By means of exact diagonalization, we analyze the evolution of ground state structures as 
the rotational frequency increases. Various kinds of ordered structures are observed. In some 
cases, hidden interference patterns exhibit themselves only in the pair correlation function; in some 
other cases explicit broken-symmetry structures appear that modulate the density. For < 10 
atoms, the standard scenario, valid for large sytems (i.e. nucleation of vortices into an Abrikosov 
lattice, melting of the lattice, and subsequent appearance of fractional quantum Hall type states 
up to the Laughlin state), is absent for small systems, and only gradually recovered as N increases. 
On the one hand, the Laughlin state in the strong rotational regime contains ordered structures 
much more similar to a Wigner crystal, or a molecule than to a fermionic quantum liquid. This 
result has some similarities to electronic systems, extensively analyzed previously. On the other 
hand, in the weak rotational regime, the possibility to obtain equilibrium states, whose density 
reveals an array of vortices, is restricted to some critical values of the rotational frequency $7. The 
vortex contribution to the total angular momentum L as a function of $7 ceases to be an increasing 
function of $7, as observed in experiments of Chevy et al. Instead, for small N, it exhibits a 
sequence of peaks showing wide minima at the values of $7, where no vortices appear. 

PACS numbers: 73.43.-f, 05.30.Jp, 03.75.Kk 
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I. INTRODUCTION 



A. Ordered structures in ultracold gases and their detection 

Ordered structures, and in particular hidden ordered structures have been a subject of 
intensive studies in the physics of Bose- Einstein condensates (EEC's) 1210,0,01, and more 
generally, in the physics of ultracold atoms. 

The paradigm example of such structures is realized in the interference of two EEC's, ob- 
served in seminal experiments of Ref. j^. Suppose that, despite the superselection rule, one 
could prepare the two condensates in coherent atomic states, characterized by fluctuating 
atom numbers A^^i, N21 but sharply defined phases 0i, 02? minimizing the Heisenberg uncer- 
tainty relation for the number and phase operators. Then, the phase difference = 0i — 02 
would determine the position of the interference fringes. Similarly, if we prepared two con- 
densates by, say, splitting a parent condensate with fixed number of atoms = A^i + A^2; we 
would arrive at sharp values of both A0 and Ni — N2- Amazingly, the interference pattern 
will also appear if the two EEC's in the Fock states (with fixed A^^i and A''2) overlap. The 
reason is, as pointed out in Ref. 0, 0], that as soon as we start detecting atoms without 
knowing which condensate they originate from, the measurement will introduce the nec- 
essary uncertainty of the atom numbers, narrowing the relative phase distribution. As a 
consequence, an interference pattern with a sharply defined A0 is obtained in each realiza- 
tion of the measurement. We may say that the measurement process uncovers the otherwise 
hidden interference pattern in the two-point first order correlation function of atomic cre- 
ation and annihilation field operators, (^"''(r)^(r')). If experimentally averaged over many 
realizations the interference pattern vanishes, since each realization leads to a different and 
completely random A0. Similar measurement induced structures, and the interplay between 
single shot and averaged results have also been discussed in the context of dark solitons in 
EEC 0. 

Other types of ordered structures occur in rotating EEC's. In the standard scenario. 



as the rotational 
structures 
Abrikosov lattice 



requency increases, more and more vortices appear in form of regular 
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. As their number grows, they organize themselves in a triangular 
. Note, that in principle the ground state of the rotating system in 



a harmonic trap should ideally be rotationally invariant and have a fixed total angular 
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momentum L, ergo it should not exhibit any structures that break rotational symmetry, as 
the Abrikosov lattice does. In reality though, the preparation of vortices is performed by a 
"laser stirring" process that breaks rotational symmetry, and introduces significant couplings 
between states with different total angular momenta [13]. Here, one deals with a situation 
in which the preparation process (which may also be regarded as a form of measurement) 
reveals elsewise hidden structures in the density of the condensate, i.e., in the one-point first 
order correlation function (\l/^(r)^(r)). 

As it is very well known from quantum optics, measurements of first order correlation 
functions (first order "coherence") do not always reveal the underlying structures. In order 
to see them, one has to measure higher order coherences, such as second order correlation 



functions (^^(ri)^^(r2)^'(r^ 
Michelson interferometry 



\l/(r4)). The paradigm example for this necessity goes back to 
14| which measures first order coherences and is sensitive to at- 
mospheric fluctuations. This deficiency of Michelson interferometry has stimulated Hanbury 
Brown and Twiss to measure the intensity-intensity correlations of the radiation coming 
from Sirius, which in turn allowed them to precisely determine the coherence length and the 
angular size of this star. 

Measurements of second order correlations play an important role in the physics of ul- 
tracold gases (for earlier works on atomic beams, see Q]). The most directly measurable 
quantity is the density- density correlation (pair correlation function, called pc function be- 
low): (^^(ri)^^(r2)^(r2)^(ri)), which formally is the two-point second order correlation 
function of the atomic field operators. This function has been directly measured in a recent 
atom counting experiment of the Orsay group 17], for the first time directly demonstrating 
an atomic Hanbury Brown-Twiss effect for thermal atoms and the second order coherence of 
a BEC. Earlier, a 4-point second order correlation function has been measured in Hannover 
[isl ] where density-density correlations of interfering condensates have been monitored in 
order to precisely determine the phase coherence length of quasi-lD condensates, in full 
analogy to the Hanbury Brown and Twiss method. 



Recently, yet another tool, i.e., noise interferometry has been proposed to analyze visible 
and hidden structures appearing in various quantum phases of ultracold gases 

HQ- This 

method also allows to determine density-density correlations, and has been used by several 
groups to study, for instance, interference of independent BEC's |2l|], residual coherence 
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22j | , and pair correlations of fermionic atoms in a Fermi 



and lattice order in Mott insulators 
superfiuid 

At this point it is necessary to mention that the double (spatial and ternporal) Fourier 



transform of the pc function is known to be a dynamical structure factor 
measurable, for instance in Bragg scattering experiments 



and is also 



B. Rapidly rotating ultracold gases 

Recently, a considerable interest has been devoted to rapidly rotating ultracold gases, 
which also exhibit various kinds of ordered structures, and should therefore be investigated 
along the lines discussed in the previous subsection. 

Typically, one considers a quasi 2D gas in the XY plane rotating around the Z axis 
with frequency Q, and confined in a harmonic trap of frequency uj±. As stated above, in 
macroscopic atomic clouds for moderate Q < uj±, the Abrikosov vortex lattice is formed 
Q. As n approaches ... the vortex lattice melts, and the system evolves through a 
sequence of strongly correlated states j26|, |2|^ . Finally, in the regime of critical rotation, it 
forms a bosonic Laughlin liquid 25^ . 

Alternatively, the various regimes of rapidly rotating gases can be described in the ter- 
minology of fractional quantum Hall effect (FQHE) theory J^]. The crucial role is played 
by the direct analog of the Landau level filling factor in the FQHE which can be related to 
the number of vortices hj u = N/N^ as defined in the BEC mean field description valid 
for large systems and moderate rotation. 

The first papers on atomic systems 3, El have considered the lowest Landau level (LLL) 
for strong enough rotation. Recently, correlated liquids at z/ = k/2 for k = 1,2,3, .. . for 
u < Ur — 6 — 10 have been discussed 



states from the Rezayi-Read 
the Moore-Read paired state 



29^ . These states resemble to a great extent the 

a'""R) hierarchy = 2 is 

etc. It has been shown that the presence of but a small 
amount of dipole-dipole interactions unambiguously makes the RR state with A; = 3 the 
ground state at filling = 3/2. This state is particularly interesting, since its excitations 
are both fractional and non-Abelian. The validity of the LLL approximation for rotating 
gases is also discussed in the recent preprint js^. 
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Most of the literature on ultracold rapidly rotating gases aims at considering relatively 
large systems and even the thermodynamic limit. In numerical simulations, either periodic 
(torus) or spherical boundary conditions are used. Unfortunately, in the N oo limit the 
gap separating the Laughlin state from its excitations vanishes. Observation of Laughlin 
states not only requires to reach the LLL, but also to control very precisely a delicate balance 
between Q and uj_i_. Despite the pro gres s in experimental studies of vortex lattices [sS. Is^. 



and first steps towards LLL physics 



35[, experiments have not yet reached this regime. 



The problems related to the short range nature of the Van der Waals forces can be over- 
come in dipolar gases, i.e. gases that interact via magnetic or electric dipole moments (for 
a review see Q]). Rotating dipolar bosonic gases are expected to exhibit exotic behaviour 
in the weakly interacting regime S^l, whereas fermionic dipolar gases have a finite gap for 
the u = 1/3 Laughlin state [.38'|. The first observation of BEG of a dipolar gas of Chromium 
atoms with large magnetic dipole has been recently reported 3^, and several groups are 
trying to realize and control an ultracold gas of heteronuclear molecules with large electric 
dipole moments |4Q]. 

Another way to create highly-correlated liquids could be, not to mimic effects of magnetic 
fields by rotation, but by appropriately designed control of tunnelling phases in optical lat- 
tices j4l|. In trapped gases, a similar effect may be realized by employing electromagnetically 



induced transparency 



42|. 



However, the most promising way towards the FQH regime and related states may be 
achieved by use of an array of rotating optical microtraps, either in an optical lattice 
or created by an array of rotating microlenses j^. 

In such arrangements, it will be natural to study mesoscopic, or even microscopic systems 
of few atoms. Such experiments demand careful theoretical studies of few atom systems us- 
ing possibly exact methods, such as exact diagonalizations of the Hamiltonian with open 
boundary conditions in the presence of the harmonic trap, or even a deformed trap. Such 
studies have recently been 

tional quantum Hall states of a few bosonic atoms have been investigated in detail. We 
continue the studies of small systems of atoms in rotating traps, and expand them in the 
present paper. 
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C. Plan of the paper 

The main focus of the present paper is to study and analyze the ground state (GS) 
ordered structures and interference patterns (IP) of two-dimensional Bose systems of few 
atoms confined in a harmonic trap and submitted to fast rotation around the perpendicular 
axis. We investigate here, on one hand, the situations where the cylindrical symmetry is 
explicitly broken so that the one particle density already exhibits ordered structures due to 
the coherent mixing of degenerated GS with different total angular momenta. Furthermore, 
we consider situations in which the ordered patterns are hidden in a pure single state with 
well defined angular momentum, and are evident only through inspection of the pair corre- 
lation (pc) function calculated by means of exact diagonalization. This formalism turns out 
to be the appropriate method to deal with small systems, for which the assumptions made 
in mean field theories do not apply. Such systems are experimentally accessible, and both, 
density and pc functions, are measurable by various experimental techniques discussed in 
the previous subsection. Experimental information on the IP can be obtained in the last 
case. 

According to our findings, the behaviour of confined systems of few atoms strongly differs 
from the behaviour of large systems. These differences are not only related to the nucleation 
of vortices in the regime of relatively slow rotation, but also to the nature of the Laughlin 
state and other highly-correlated states, when the rotational frequency is close to the trap 
frequency. 

In particular, we obtain that for < 10 atoms the standard scenario valid for large 
sytems (i.e., the nucleation of vortices into an Abrikosov lattice, melting of the lattice, and 
subsequent appearance of fractional quantum Hall type states up to the Lauhglin state), 
is absent for small systems, and is only gradually recovered as A^ increases. On the one 
hand, the Laughlin state in the strong rotational regime contains ordered structures much 
more similar to a Wigner crystal or a molecule than to a Fermi liquid. This result has some 
similarities to electronic systems, extensively analyzed previously. On the other hand, in the 
weak rotational regime, the possibility to obtain equilibrium states whose density reveals 
an array of vortices is restricted to some critical values of the rotation frequency Q. The 
vortex contribution to the total angular momentum L as a function of Q ceases to be an 
increasing function of fl, as observed in experiments of Chevy et al. pLj. Instead, for small 
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A^, it exhibits a sequence of peaks showing wide minima at the values of where no vortices 
appear. 

This paper is organized as follows. In Sec. II we describe our system (Section IIA) and 
address the questions related to the realization and analysis of ground states with hidden 
(Section IIB), or explicit (Section IIC) broken cylindrical symmetry. In Sec. Ill, the main 
results of this work are presented. Finally, in Sec. IV, we compare our findings with previous 
results in the literature, and draw our conclusions. 



II. ORDERED STRUCTURES IN GROUND STATES: BROKEN CYLINDRICAL 
SYMMETRY 



A. Description of the system 

Our system consists of bosonic atoms trapped in a rotating parabolic potential. The 
Hamiltonian in the rotating reference frame reads 



N 



i=l 



g^S{fi-fj), (1) 



i<j 



where r = (x, y), uj± is the trap frequency. A* = x r is the vector potential, z is the 
unitary vector along the Z direction and B* = V x A* = "^M^^ jg ^j-^e effective magnetic 
field of an equivalent system of electrons submitted to a magnetic field perpendicular to the 
XY plane (we use here the symmetric gauge). Thus, the rotation of the trap has formally 
the same effect on atoms of mass M , as a magnetic field has on electrons; the electronic 
charge — e and the speed of light c are solely introduced for reasons of algebraic equivalence. 
V = g Yli<j is the contact interaction potential, where g is the interaction coefficient 

that approximates the potential of the Van der Waals forces between the atoms in the very 
dilute limit. We assume the rotational frequency to be large enough to restrict the system 
to the lowest Landau level (LLL) regime, and choose the appropriate Fock-Darwin single 
particle (sp) wave functions with no nodes in the radial direction, as the basis in order to 

n 

represent all operators |47l |. 

I m) = ^= {\Y e"l^l^/2^^ (2) 
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with ^ — y 2Mu)j_ ' generalized complex coordinates z — x + iy. 
The Hamiltonian can be written in second quantized form as, 

H = aL + pN + V, (3) 

where a = h{u)± — fl), /3 = huj±, L and N are the total z-component angular momentum 
and particle number operators, respectively, and 

V = ^ ^ ^1234 a[ala4 03 , (4) 

mim2m3m4 

where the matrix elements of the interaction term are given by 

/ I TA I \ 9 5mi+m2,m3+m4 [mi + 7712)1 , . 

Vi2u = mi m2 y m3m4 = = . 5 

Here, the operators a| and Oj create and annihilate a boson with single-particle (sp) angu- 
lar momentum m^, respectively. The cylindrical symmetry of the Hamiltonian allows the 
diagonalization to be performed in different subspaces of well defined total z-component of 
angular momentum L = ^^i mj. 

Fig. 1 shows the total angular momentum of the GS of a system oi N — 5 particles 
while Q grows from zero to u!±, the maximum possible value before the system becomes 
centrifugally unstable. We observe that the GS angular momentum remains constant for 
a finite range of Q until transitions to new angular momenta take place at critical values 
labelled as flcn- Not all L-values can be associated with the GS's. However, on the steps, 
different L-states may be degenerate in energy as in the case of the states with L = 8, 10 and 
12 on the third step (indicated by stars in Fig. 1). The last possible GS at L = N{N — 1) 
is the Laughlin state, for which the interaction energy is zero due to the fact that the wave 
function of each atom has zeros of order two at the positions of the other N — 1 atoms; this 
can be easily deduced from the analytical expression of the many-body wave function given 

by, 

^LaugkUn = J^U^'i - ^.^^e" ^ l^'l ^^^^ . (6) 
i<j 

In the following, the GS of the system in the L-subspace is denoted by ^l- In order to 
study the nature of the GS's, it is useful to analyze the expectation values of some relevant 
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FIG. 1: Change of the GS angular momentum Lgg for = 5 as the rotation frequency increases; 
transitions take place at critical values of the rotational frequency labelled by f^cn- (a = h{uj± — Q)). 

operators. First of all, it is crucial to realize that the density operator defined in first 
quantization as 

N 

p{r) = ^S{r-ri) (7) 

i=l 

does not exhibit any interference pattern when calculated for a definite "^l, as it can be 
inferred from its analytical expression in second quantized form 

Pi^ = J2^Hr') I Sir- P) I (PjiP))ala„ (8) 



where \4>i{r)) = |mj) as in Eq. Due to angular momentum conservation, the operator 



ajaj selects only one sp state and, as a consequence, it loses all information contained in 
products of different amplitudes, thus losing the interference pattern. It solely preserves the 
information of individual densities, e.g., 

N 

p{r) = {^L I P(r) \^l)=J2\ I' ' (9) 
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where Oq is the total occupation of the sp state |mj) in the GS. In effect, p{r) can only 
represent cylindrically symmetric distributions. However, we note that this cylindrical sym- 
metry is a direct consequence of the definition of the operator p(r), and is not necessarily a 
manifestation of the symmetric nature of the GS. 

To exhibit ordered patterns and analyze the GS structures, we proceed in two different 
ways; one investigates the pair correlation function for states with fixed L (Section IIB), the 
other combines different \1/l's (Section IIC). 



B. Ordered structures in pair correlation functions 

In order to analyze the internal structure of relevant states, we consider the following 
operator 

N 

pir,ro) = ^6{fi- fo)6{fj - f), (10) 
i<j 

which yields the conditional probability to find an atom at r*, when another is simultaneously 
found at tq. This operator contains information that originates from the amplitudes of sp 
wave functions, and not only from their density as it has been in case of the single particle 
density operator. In second quatized formalism, its expected value with respect to^i reads 

p{r,ro) = j^a*ap'0*(f)(/>*(fo)(/>fe(r)0i(fo)($p | ala]aiak \ (11) 

ijkl pp' 

where 

^L = J2^p%, (12) 

p=i 

and $p are the bosonic Fock A^-body states of the basis in the L-subspace of dimension 71^. 
The condition i+ j = k + l must be fulfilled for reasons of angular momentum conservation. 
It should be stressed that p(r, ro) in Eq. |Tn|l obviously differs from the single particle density 
matrix 

nW(f,r') = (^+(f)^(r"^)), (13) 

which defines the off-diagonal long-range order that characterizes Bose condensation j^. 
The operator p is a two-particle operator, whereas n'^^^(r', r') = ^"'"(r)^(r') is a single- 
particle operator; in particular, p(r) = n^^\r,r), whereas p(r) = J dfop{r,ro). As a 
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rule of thumb, if n^^\r,r') reveals symmetry breaking, so does p, whereas the opposite is 
not necessarily true. 

Eq. (fTT|) can be interpreted as the sum of products of amplitudes at r weighted by a 
factor that depends on rg, and on the GS via the ap coefficients. In the particular case 
•tq = 0, cylindrical symmetry is recovered, since in this case / = j = is the unique non-zero 
contribution. This implies i = k and yields 

p(f,o) =1 0o(o) p Yl I I' E^p"^''^-)' (14) 



pp 



which is independent from 6. 

In order to understand the role of the parameter fo in | /3(r, ro) | ^I/l) as a function 
of f, we consider 

1 ^rrin mj mi 

TT^ V'^n! ^/mjl y/mk\ y/m' 

in units of A. As I— j = n — k follows from angular momentum conservation, the angular 
dependence reads 

Evidently, if ro is fixed, the change of 9q is nothing but a rigid rotation of the function. 
In other words, any arbitrary choice of Oq fixes the origin of angles, and breaks cylindrical 
symmetry, in the analogous way as it happens in experiments which perform a single shot 
measurement. Within this point of view, the experimental measurement and the choice of 
Oq are equivalent processes (see for instance [30|)- 

The expected values of the pc function for tq 7^ can reveal very different situations: 
from circular symmetric structures showing no spatial correlation, to ordered structures that 
reveal intrinsic Wigner molecules or crystals while passing through all possible intermediate 
states, as it is shown in Section III. 



C. Ordered structures in the density: superpositions of different L-subspaces 



In this subsection, we consider ordered structures in GS's with no well defined angular 
momentum in two different situations. On the one hand, we build (somewhat ad hoc) linear 
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combinations of different \I/l's to explicitly reveal the structure present in the expectation 
value of the density operator. We obtain ordered IP's for combinations, whenever one of 
the contributing \E'l's has an ordered hidden IP contained in its pc function. On the other 
hand, after the introduction of an anisotropic term to the Hamiltonian which mimics the 
deformation introduced by the stirring laser, we perform numerical diagonalization without 
the restriction of angular momentum conservation, and obtain in this way GS structures 
with broken symmetry. These exact calculations give hints, how to construct approximated 
GS superpositions in the previous ad hoc construction. 

To explain more precisely, how the first procedure works, we start from the pc function 
and observe what kind of ordered structures can be expected. In case of the Laughlin state 
L = 20, = 5, the pc function suggests that the atoms form a pentagon. Quite generally, 
the best way to visualise this structure within the first procedure is to form a superposition 
A^'x+i?\l'i_|_jv, where is the GS that contains a hidden ordered IP. It is easy to understand 
this result from what follows. The terms that contribute to the broken cylindrical symmetry 
are those of the form {"^l \ o-l'^j I ^l+m) with M > 1. However, it is necessary to arrive at 
M = in order to obtain contributions from all the sp states contained in "^l- To be more 
precise, none of the combinations from L = 20 + 21 to L = 20 + 24 reproduces the structure 
of L = 20 for = 5. It is necessary to combine L = 20 + 25 to obtain the regular pentagon 
implicit in \E'2o- The best contrast is obtained for A = B = 1, and not for a small amount 
of "ifL+N as one would expect if only a perturbation would be necessary. 

It is important to stress, that the ordered hidden IP was obtained in the Laughlin state 
L = 20 for = 5. In order to assure that this state is the GS, a very small amount of kinetic 
energy is necessary in such a way that the states L + M for M > are quasi-degenerate. 
Then, the combination considered above corresponds to a "legitimate" GS. 

The second procedure, followed in order to obtain ordered structures as, e.g., multiple 
vortex states, was suggested by experiments. In the experimental setup described by Chevy 



condensate rotating under the action of a stirring laser that produces anisotropy in the XY 
plane. Subject to this anisotropic potential, the state with vortices is a GS, and it survives 
during the time of fiight (TOF) detection as an excited state of the restored symmetric 
Hamiltonian after the trap is switched off. 




vortices are generated as equilibrium states of a Bose 
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With this idea in mind, we introduce an additional anisotropic term in the Hamiltonian 
given hj Vp = A J^iLii^"^ ~ y'^)i or i^i second quantized form as j3l 

Vp = ^ jA/m(m - l)a+a„„2 + ^/(m + l)(m + 2)a]„a„+2 • (17) 

m 

We assume this term to be a small perturbation of the system, thus, ^^'^'^^^^ << 1, and 
perform exact diagonalization to obtain the GS of H + Vp. 

Amazingly, the structure of, say, two vortices can only be obtained at very specific plateau 
steps Qcn- There, the GS is a combination of quasi-degenerate \l/2,-states which are coupled 
by a perturbing term slightly larger than their energy difference but much smaller than 
the next eigenenergy. As a consequence, the linear combination of the states above has 
nearly equal coefficients. Thus, a direct ad hoc combination of degenerated states of the 
symmetric Hamiltonian of Eq. are educated guesses to reveal underlying structures. 
This combination was previously used by Wilkin et al. in exact diagonalization calculations 
to obtain two vortices j26|. More precisely, the unique situation where vortices are generated 
in the density corresponds to the steps in the Lgg dependence on fi, where a degeneracy of 
states with different L takes place at Vlcn- At first sight, this result does not agree with 
the experimental results reported by Chevy et al. However, it can be attributed to 

a essentially different behaviour of systems with a large and a small number of atoms, 
respectively. As grows, the size of some of the plateaus shown in Fig. 1 drastically 
shrinks in such a way that finite ranges of fi-values with energetically degenerate states 
become possible; not only at critical values ficn- In Fig- 2, we show the appearance of such 
microplateaus obtained for A^ = 6, 7, 8, and 9. 

Deduced from these features, our prediction is that the experimental graph, analogous to 
the one displayed by Chevy et al. |^ in their Fig. 2 which shows a monotonous growth of the 
vortex contribution to Lgs as a function of fi, from the first vortex nucleation (at fici) to the 
turbulent regime, would look radically different for small A^. We expect that it would present 
a curve with minima at those values of where the GS is deep in a plateau, and exhibits 
no vortices in the density. Vortices will solely be visible on microplateaus surrounding VL^n- 
For fixed but larger A^, the microplateaus contain more and more states as Vl increases, and 
thus a larger number of vortices is nucleated. Ultimately, in the large A^ limit, the number 
of vortices becomes proportional to the rotational frequency. 
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Number of particles/ Angular momentum of the ground state 

* N=9/0,9, 14, 16,21 ,24,28,33,35,37,44,56,72 

^ N=8/0,8, 1 2, 1 4, 1 8,24,30,32,35,42,56 

▼ N=7/0,7,1 2,1 5,1 8,24,30,35,42 

^ N=6/0,6,1 0,1 2,1 5,20,24,30 

■ N=5/0,5,8,1 2,1 5,20 

▲ N=4/0,4,8,12 

N=3/0,3,6 



I 

0.35 0.3 0.25 0.2 0.15 0.1 0.05 

n (CO^-Q) 

FIG. 2: The same as in Fig. 1 for = 3 to 9 from bottom to top. Graphs are vertically shifted 
for clarity. The (quasi-)degenerate states within the same step are not included. 
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III. NUMERICAL RESULTS 



In what follows, we display the results obtained from exact diagonalization for = 1 in 
units of hujj^. The values uj_[_ = = were considered for density and correlation functions 
of single states, or combinations of states with well defined angular momentum, as the 
diagonalization of H depends only on the interaction energy. In contrast, specific values 
of a and (3 (see Eq. (jS))) are considered when the diagonalization of if + is performed. 
Correlation functions are always displayed in pairs, a 3D and a contour plot, unless otherwise 
specified. We consider A and huj±_ as units of length and energy unless otherwise stated. 

For ^ 1, the first vortex is nucleated at fici = ^± ~ when the transition from L = 
to L = N takes place at the first step of Lgs{Q). The L = N state whith a vortex at the 
centre is the GS from Qci to Qc2 (see Fig. 1). This state is characterized by a high occupancy 
of the m = 1 sp state and circular symmetry possessing no space correlations. In contrast, 
for small values of a vortex is not clearly manifested unless a considerably large number 
of atoms is considered, as it is shown in Fig. 3, where the pc of the L = N state is displayed 
for = 3, 10 and 20. In all the cases, the density has a minimum at the origin and tq is set 
to its maximum. The slow tendency to recover the behaviour of large condensed systems is 
explained in Fig. 4, where the occupation of the m = 1 sp state over is shown. Namely, 
for few atoms, the angular momentum of the state is not fully due to vortices. For a large 
number of atoms, Nakajima and Ueda have analyzed the formation of the first vortex, from 
its initial nucleation at the cloud boundary towards its final stabilization at the centre 
our results share some similarity with this situation: for small A^, the vortex is not yet fully 
inside the trap, and as A^ increases, it approaches the trap centre from the boundary of the 
cloud. 

In Fig.'s 5 to 9, we show the main results for A^ = 3 including the ground state evolution 
as Q increases. Fig. 5a shows the lowest eigenenergies for each L, the so-called Yrast line. 
The initial points of the plateaus, at L = 0, 3 and 6 are the unique possible GS's (besides 
degeneracies at the steps). A general result is that the plateau previous to the Lauhglin 
state (from L = 3 to L = 5 in this case) has always A^ points. Fig. 5b is similar to Fig. 1 
for A^ = 5. The densities from L = to 9 are shown in Fig. 6. The parameter tq used 
in the pc calculation was set to the maximum of the density unless it is located at the 
centre. Then, tq = 1 is used if not stated elsewise. Fig.'s 7 and 8 display the pc function 
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for L = 0,3,4,5,6, and 9. The system evolves from a completely "condensed" system at 
L = to the Laughlin state at L = 6 where a clear triangular structure appears. The loss 
of condensation is related to the increase of space correlations. To complete the analysis, 
we show in Fig. 9 the evolution of sp occupations, and demonstrate that "macroscopic" 
occupation of a specific sp wave function vanishes as L increases. 

In order to see how the previous general tendency evolves as N increases, we analyzed the 
= 5 case. In Fig. 10 we show the densities of the GS's from L = to the Laughlin state, 
and in Fig.'s 11 and 12 we display their pc. The same tendency towards space ordering 
at the Laughlin state is clear. In addition, it can be inferred from the L = 20 case that 
correlations are stronger for nearest neighbours as a manifestation of partial long range order 
in finite systems. In Fig. 13 the occupations of sp states are shown. It is remarkable that 
some indications of the Laughlin state typical for large systems, i.e., a fiat density at the 
central part and a hump at the edge are already manifested in such small systems, as it can 
be seen in the last graphs of Fig. 13 and Fig. 10, respectively. Moreover, the density at the 
origin is very close to l/(27r), as necessary for a homogeneous system at filling factor 1/2. 

In Fig.'s 14 and 15, we concentrate on the Laughlin state for = 3, 4, 5, 6, 7, and 8 atoms. 
The left hand side picture for = 3, 4 and 5 contains the pc function and the right hand 
side displays the density of the superposition made from L and L + N. This superposition 
of quasi-degenerated states becomes a possible realization of a GS for a certain value of Q, 
as commented previously. The value of the parameter vq in the pc functions was obtained in 
a different way as the one used previously. Taking advantage of the fact that the Laughlin 
wave function (Eq. is the exact solution, and knowing that its pc shows a ring shape 
structure of an unknown radious vq, we can maximize the probalibity distribution given by, 

I "^Laughhnirl, r^, ...r^) p= e"'^ (18) 

where 

i i<j 

(with q = 2 for the bosonic Laughlin state), or equivalently, minimize T with respect to tq. 
Minimization yields tq = ^/W^^ (or tq = if one atom is at the origin as for = 6, 7 
and 8) which is always smaller than the size of the system given by i? = yAN^^2. To see 
the evolution of these hidden ordered structures as A^ increases, the "degree of correlation" 
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C as a function of is displayed in Fig. 16; C is defined as the height of the maximum 
peak in the pc function. It decreases for increasing as it is expected in order to recover 
the quantum hquid character of the Laughhn state for large systems. 

Finally, for = 5 and 6, Fig. 17 shows patterns of two incipient vortices obtained from 
full diagonalization oi H + Vp at h{uj±_ — Vt) = 0.0458 and h{uj±_ — Vt) = 0.05904 in units of 
huj^, respectively. The possibility to obtain these patterns for the given anisotropy strongly 
depends on the possibility to obtain truly degenerate states with angular momenta L and 
L ± 2 at one of the steps Lgs{Q). For N = 7 we were not able to find such steps. For = 5 
the result of the diagonalization at the third step, where L = 8, 10 and 12 are involved, is 
that the weights of L = 8 and 10 within the expansion of the GS are much larger that the 
weight of L = 12. This leads in effect to a state with expected vortex angular momentum 
lower than 10, in agreement with the results demonstrated in Ref. 50| related to the fact 
that the contribution of a vortex to the total angular momentum depends on its distance 
from the origin, it runs from A^ at the centre to zero at the trap boundary. 



IV. DISCUSSION AND SUMMARY 

The main result of this paper is the identification of important differences between large 
and small systems of rapidly rotating cold bosonic atoms. These differences can be under- 
stood by looking at the expected values of the density and the pc functions, on which we 
have concentrated our analysis. The characterization of small samples in rotating traps has 
recently attracted increasing interest due to the possibility to deal with a few number of 
atoms per well in optical lattices 45^ . 

Within the regime of low rotational frequency, we obtain that a relatively large number 
of atoms is necessary to nucleate the first vortex carrying A^ units of angular momentum. 
The evolution towards the condensed state with L = A^ is shown by the increase of the 
occupation of the sp m = 1 state as A^ increases. On the other hand, in the regime of strong 
rotation, space correlations increase significantly, and in the Laughlin state a hidden ordered 
structure modulates the pc pattern. For small systems, the atoms sit around a ring of radius 



To = — 1 [N = 3, 4 and 5) or tq = vN {N = 6, 7 and 8). The degree of correlation 
defined as the height of the peaks in the ordered pattern decreases with A^, evolving towards 
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a non-correlated structure of a quantum liquid. We have argued about the observability of 
the ordered IP in similar experiments as those reported, for instance, by Foiling et al. j2^. 

N^umerous references have analyzed the Wigner structures of few electrons 
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mostly for filling factors less than 1/3, and exhibited in the density (and not in the 
pc function). It is important to remark that it is well established from exact diagonalization 
studies in a torus geometry and from the analysis of the Laughlin wave function, that states 
of filling factor 1/3 for electrons and 1/2 for bosons are fermionic quantum liquids in the 
thermodynamic limit (cf. Whenever the Laughlin function is a good approximation, 

its implied properties are independent from the interaction. However, for small confined 
systems the previous results do not apply, and the analysis of some of their properties relies 
on the competition between the kinetic and interaction energies, aside from their statistics. 

Suggested by our results, a possible explanation for the realization of Wigner molecules in 
the Laughlin states for few bosons is the following. The first observation is that the nature 
of the GS does not depend on the kinetic part as the diagonalization is fully determined by 
the interaction, in other words, the structure does not result from the competition between 
different kinds of energy. In addition, as the repulsive interaction energy is zero in the 
Laughlin state, it seems that the reason why the atoms choose symmetric and well separated 
positions is due to two conditions, firstly, the system must have a large angular momentum 
given by L = A^(A^ — 1) (which means large distances from the origin) and secondly, each 
atom is surrounded by a quasihole (which leads to effective mutual repulsion). This last 
statement is supported by the following observation. The contour plots of = 3 for L = 4 
and L = 5 in Fig.'s 7 and 8 suggest that in those precursory states (the Laughlin state 
has L = 6), quasi- holes not attached to atoms are created without cost of internal energy, 
the contribution to the angular momentum of each one would evolve as 1/3 (in L = 4), 
2/3 (in L = 5) and 3/3 until the Laughlin structure becomes possible with one quasi-hole 
attached to each atom, lowering the interaction energy. A final observation relates to the 
evolution of this behaviour as increases. Due to the fact that the dependence of L on A^ 
is quadratic, the increase of L with A^ is more efficient for large A^, and atoms do not have 
to be widely separated. Thus, the symmetric distribution tends to disappear. Recently, 
it has been proposed that the phenomenology of strongly correlated bosonic and fermionic 
rotating systems converges to te case of classical particles, and finally cristallizes at high 
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rotational frequencies j59|. Our results exhibit some traces of such cristallization, but it 
should be pointed out that this "cristallization" ceases to be manifested when the bulk 
structure starts to dominate the system. This happens for bigger particle numbers, where 
the GS at Laughlin angular momentum will start to behave more and more like a true 
quantum liquid as pointed out above. 

Finally, precursors of two- vortex arrays are obtained as the ground states of an asymmet- 
ric Hamiltonian that models the experimental setup used to increase the angular momentum 
of a trapped Bose condensate by a stirring laser. We conclude that the possibility to nucle- 
ate vortex patters in the density is restricted for small to the specific values Qcn at the 
steps where several degenerate or quasi-degenerate states of different angular momentum 
L coexist. This produces peaks in the vortex angular momentum dependence on Q. We 
predict these peaks to broaden as increases, due to the appearance of "microplateaus" 
which in turn lead to finite ranges of Q values where quasi-degenerate states coexist. 
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FIG. 3: Pair correlation function for A?^ = 3, 10 and 20 of the L = N state. The parameter rO is 
equal to 0.8, 0.95 and 1.0 in units of A respectively. 
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FIG. 5: For N = 3, a) interaction energy as a function of total angular momentum (Yrast line), 
b) angular momentum of the GS over a. The critical values for a at the steps are: 0.1194 and 
0.0398 in units of nuj±. 
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FIG. 6: N=3, density of the L-states 
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FIG. 8: The same as Fig.7 for L = 5, 6 and 9, rO = 1.0, 1.3 and 1.7 respectively. 
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FIG. 9: For A'^ = 3, total occupations of the single particle states of angular momentum m for 
several L-states. 




FIG. 10: N=5 density of the L-states (GS). 
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FIG. 11: Pair correlation of the L-states (GS) for AT = 5 and L = 0, 5, and 8. rO = 1.0, 0.9 and 
1.0 respectively. 
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FIG. 13: N=5 occupations of the L-states (GS). 
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FIG. 14: For N = 3, 4 and 5 the left hand plots show the pair correlations (contour plots) for 

L = 6, 12 and 20 respectively (the Laughlin states), with rO = \/N — 1. The right hand 3D-plots 
show the density of the mixtures (with equal coefficients) of angular momenta 6 + 9, 12 + 16 and 
20 + 25 respectively. 
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FIG. 15: For N = 6, 7 and 8 pair correlations ( 3D-plot and contour-plot) of L = 30, 42 and 56 

respectively (the Lauglilin states). We consider rO = y/N. 
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FIG. 16: Correlation degree (see text) of the Laughlin state as a function of N. 
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FIG. 17: For N = 5 and 6 density ( 3D-plot and contour-plot) of the two vortex structures. For 
iV = 5, L = 8 + 10 + 12 and for iV = 6, L = 10 + 12 + 14. 
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